Optical recombination of biexcitons in semiconductors 
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We calculate the photoluminescence spectrum and lifetime of a biexciton in a semiconductor 
using Fermi's golden rule. Our biexciton wavefunction is obtained using a Quantum Monte Carlo 
calculation. For a recombination process where one of the excitons within the biexciton annihilates, 
we find that the surviving exciton is most likely to populate the ground state. We also investigate 
how the confinement of excitons in a quantum dot would modify the lifetime in the limit of a large 
quantum dot where confinement principally affects the centre of mass wavefunction. The lifetimes 
we obtain are in reasonable agreement with experimental values. Our calculation can be used as a 
benchmark for comparison with approximate methods. 



I. INTRODUCTION 

Information about many-body systems of electrons and 
holes is primarily obtained through their luminescence. 
While good calculations on the ground states of these 
systems have been performed for a variety of electron and 
hole densities, good estimates of luminescence require the 
inclusion of final excited states, for which calculations are 
difficult. Thus, it is important to quantify how excited 
final states affect luminescence. In this paper, we solve 
essentially exactly the problem of decay of the simplest 
correlated system, namely the biexciton. 

Luminescence of these biexcitons, which are bound 
states of excitons (electron- hole pairs), is potentially 
relevant for studying luminescence from semiconductor 
lasers. While semiconductor lasers correspond to the 
dense limit of the electron-hole phase diagram, where 
the electrons and holes ionise into a plasmai, it is well 
known that the inclusion of correlation effects with the 
surrounding electrons and holes is necessary for a good 
estimate of semiconductor luminescence^Ti^. Moreover, a 
variational Monte Carlo study of the electron-hole phase 
diagram has revealed a surprisingly large excitonic insu- 
lator phase, which extends well into these high densities'^. 
Thus, luminescence from excitonic states also contributes 
to light emission in the plasma. However, it is not a pri- 
ori obvious whether emission from an independent exci- 
ton can adequately describe emission from a plasmonic 
state where interaction effects are crucial. This ques- 
tion of interaction can be addressed by considering the 
biexciton as the most complicated correlated state that 
can still be treated exactly. Since the detailed form of 
the bound state wavefunction can significantly affect the 
overlap with final states it is important to go beyond the 
use of simple analytical variational wavefunctions for the 
biexciton. Instead we use a biexcitonic wavefunction ob- 
tained from a quantum Monte Carlo calculation. Thus, 
our results can be used as a benchmark for larger and 
more complicated systems that currently still require ap- 
proximate methods. 



The concrete aim of this work is twofold: we primar- 
ily want to calculate the emission spectrum of the biex- 
citon, using an essentially numerically exact biexciton 
wavefunction. We then want to find out how well the 
lifetime we calculate from this spectrum compares to ex- 
perimental and other theoretical lifetimes. 

For the calculation of the lifetime we focus on radia- 
tive recombination involving the emission of a photon. 
This process could potentially excite the surviving exci- 
ton into an excited state. In order to obtain both the 
emission spectrum and the lifetime we follow the proce- 
dure by EUiotli^ and use Fermi's golden rule, which re- 
quires the calculation of matrix elements between these 
exciton excited state and the full biexciton wavefunction. 
The wavefunctions and matrix elements are clearly dif- 
ferent for different mass ratios. Hence, we consider dif- 
ferent mass ratios where the hole mass rrih is equal to or 
much greater than the electron mass rrie for the spectra 
and several intermediate mass ratios for the calculation 
of the lifetime. 

It has been shown in previous calculations^ that the 
inclusion of the accurate semiconductor eigenstates, i.e. 
polariton states, is important for an accurate estimate 
of the biexciton lifetime. The importance of polariton 
effects for the biexciton decay comes from the fact that 
since most binding energies of excitonic molecules Sm are 
much smaller than the polariton splitting parameter fic, 
the dispersion of the final states after the decay process is 
strongly modified compared to the dispersion of the free 
exciton and photon states. Hence we estimate lifetime 
including polariton effects; in doing this we use insights 
from our calculation in the absence of polariton effects in 
order to simplify the calculation in presence of polaritons. 
We present a general formula of the typical non-radiative 
biexciton lifetime in the semiconductor, which should in 
principle work for materials where the polariton effect is 
weak. 

To the extent that our Monte Carlo biexciton wave- 
function is a good wavefunction, the remainder of our cal- 
culation is essentially exact for system where the polari- 
ton effect is negligible. We expect that shape of our emis- 
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sion spectra, calculated without polariton effects, will not 
be affected drastically even for materials where the effect 
is important. Thus, our calculation can be used as a com- 
parator for approximate methods, which could, however, 
reach non-zero densities. 

While our spectra are for 3D bulk only, we extend our 
lifetime calculation to a 3D quantum dot, since more ex- 
perimental data is available there. We compare our life- 
time estimates for a range of materials to experimental 
and other theoretical estimates. A similar procedure for 
the calculation of lifetimes as presented in this paper may 
also be applied to quantum well systems. 

In this paper we calculate the simplest optical recom- 
bination process. The question of whether this is in fact 
responsible for the exciton lifetime is however more com- 
plicated as many other mechanisms can play a role. In 
semiconductor bulk systems these mechanisms for exam- 
ple include recombination involving impurities, or traps'^, 
or Auger recombination^"^^. We discuss some of these 
later; however as our lifetimes compare reasonably with 
experimental values it seems that radiative recombina- 
tion is the main channel for biexciton annihilation. 

The paper is organised as follows: Section II discusses 
photoluminescence from biexcitons in a semiconductor, 
and section III discusses their lifetime. In section Hi Al we 
introduce the approach used to determine the biexciton 
wavefunctions and how to calculate photoluminescence 
spectrum, which we present and discuss in section IIIBI 
Section IIII Al presents the calculation of the biexcitonic 
lifetime with polariton effects and optionally confinement 
and section IIII Bl compares biexciton lifetimes for five 
semiconducting systems to experimental and some other 
theoretical values, and discusses their trends with elec- 
tron/hole mass ratio and confinement. We conclude in 
section HVl 



II. PHOTOLUMINESCENCE 
A. Methods 

In this section we present the calculation of lumines- 
cence from the biexciton. We first discuss the method- 
ology of this calculation, starting from the light- matter 
interaction Hamiltonian, and the form of the biexcitonic 
wavefunction, before quoting the general formulae for the 
relative transition rates into different excited excitonic 
states. We then present the results and discuss our ob- 
tained emission spectra. 

The light-matter interaction in second quantisation is 
given by 
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where pcv is the momentum matrix element for a 
transition of an electron with charge e and mass TOq 



between the valen ce and the conduction band and 
Aq = y/h/ {2eujkV) is the standard vector potential field 
strength. This expression contains uJk as the frequency of 
the emitted photon, the dielectric constant of the mate- 
rial e, and a unit of box quantisation volume V which can- 
cels in the latter part of the calculation. The operators 
ct „ create holes and electrons with momentum kh and 
ke and with spin quantum number a = ±f , ^, respec- 
tively, ajj. ^ creates a photon with polarisation cr' = ±1. 
Pcv is the optical or momentum matrix element between 
valence and conduction band. The factor 7 in equation 
(H)) accounts for the overlap of electron and hole spin 
states if spin-orbit coupling is neglected then this be- 
comes simply a factor of twc^^. We use the values from 
Ref. [H. 

The luminescence is determined by decay from all 
possible bound states in the semiconductors into other 
bound states through the emission of a photon. The 
rate of generating a photon per unit photon energy corre- 
sponds to the expectation value of the number of photons 
after a time t: 



k,cr 



t 



Sihoj^ E{nu,a)), (2) 



where uj is the frequency corresponding the the emitted 
photon with wavevector k, and E{hu!, a) is the dispersion 
of the photon in the medium. In this article, we consider 
luminescence from a biexciton in its ground state as a first 
approximation to understand luminescence from bound 
states in semiconductors; the result will hopefully allow 
us to estimate to what extent bound states will need to 
be considered, or whether the consideration of individual 
excitons suffices. 

In order to obtain a more explicit equation for the lumi- 
nescence than equation ([2]) , we first of all need to discuss 
the wavefunction of our biexcitonic ground state with re- 
spect to which the expectation value of equation (0) is 
calculated. Previous calculations of luminescence from 
biexcitons have used estimates or simplified versions for 
the biexciton wavefunction^i; In contrast, the biexciton 
wave function ^'be we use is of the form 



*be(R) =exp [J(R)]^'s(R) , 



(3) 



where J(R) is a Jastrow factor of the Drummond-Towler- 
Needs form-"^ containing only two-particle correlations, 
and 

\E'S(R) = (/>l(ri3)0i(r24)'/)2(?'l4)02(?'23) 

where r^j is the distance between particles i and j, with 
electrons being particles 1 and 2, and holes being parti- 
cles 3 and 4. The pairing orbitals (pnir) are of the form 



(knir) = exp 



Pl,n{p2,7i + r) P2M+r 



(4) 
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where pi^n and p2,n are optimisable parameters and F 
is a constant such that the Kato cusp conditions'^^ at 
electron-hole coalescence points is reproduced by 
electron-electron and hole-hole cusps are introduced via 
the Jastrow factor. 

Given a trial wave function 5*, assumed to be real for 
simplicity, the variational Monte Carlo (VMC) method 
is capable of evaluating the variational estimate -Eif of 
the ground-state energy Eg, 



/vl/(R)g(R)vl/(R)dR ^ 
/|v&(R)|^dR - 

by evaluating the local energy, 



^l(R) = 



ij'(R)^(R) 

*(R) 



(5) 
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at a set of M configurations {Ri} distributed according 
to the square of the trial wave function |^'(R)|^, and 
evaluating the average 



1 

1=1 



(7) 



where the "approximately equal" sign refers to the statis- 
tical error due to the finite size of the sample. In addition, 
the VMC method optimises the parameters in the trial 
wave function by means of minimizing E\ with respect 
to the parameters, using techniques such as the modified 
linear least-squares method developed by Umrigar^-'^. 

The above trial wave function is then optimised within 
VMC for different mass ratios using the CASINO code2^. 
The VMC energies are presented in Table HI where we 
also report the variance of the local energies cry, which 
is an additional measure of the quality of a wave func- 
tion. The variational energies obtained at all mass ratios 
contain a large fraction of the exact binding energy, indi- 
cating an accurate description of these systems. A com- 
parison of VMC wavefunctions with both experimental 
and theoretical values was performed in^^ for quantum 
well structures. 
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1837.36222 -1.1502(4) 0.131(1) 

183.73622 -1.1335(2) 0.01831(5) 

18.37362 -1.0338(1) 0.01587(5) 

1.83736 -0.66725(8) 0.00511(2) 

TABLE I. VMC energies and variances of the local energy for 
different hole-to-electron mass ratios. 



Using the biexcitonic wavefunction calculation with 
VMC as our ground state, we can now transform equa- 
tion ^ into Fermi's golden rule for the transition 
rate W{fujj) per energy range of a particular photon~. 
Fermi's golden rule corresponds to a perturbative treat- 
ment of the interaction^^, and has already been used by 



ElliotSi for a calculation of exciton luminescence. The 
transition rate per energy range requires a sum over all 
final states of both exciton and photon. 



k,K,{rii} 



l(*E,{n.},K|i?'l*BE)|' (8) 

(5(6; — e/(k, {rii}) — haj)S{huj — ?ic|k|). 



As final exciton momenta K are constrained by mo- 
mentum conservation, the corresponding sum disappears, 
and the momentum of the surviving exciton is — k. The 
final excitonic state with wavefunction ^'E,{ni},K is la- 
belled by the centre of mass K and the various eigenstates 
of the exciton {«;}, which are standard hydrogenic solu- 
tions. The quantum numbers {rii} correspond to n, I, 
and m for the bound and k, I, and m for the continuum 
states. As the biexciton is in its ground state, we obtain 
from angular momentum conservation that I = m — 
for the surviving exciton. is the ground state biexci- 
ton energy, and e/(k, n) the energy of the final state. 

We now insert H' from equation ^ into equation ([5]) 
and notice that it is necessary to perform several real 
space integrals in order to simplify the matrix element. 
After a separation of centre of mass and relative coor- 
dinates, the relative exciton wavefunction (f)E.n can be 
expressed in terms of a single coordinate ri = r^i — r^i, 
where Vei and r^i refer to the coordinates of electron 
and hole of the surviving exciton. The relative biex- 
citon wavefunction ^be depends on this coordinate ri, 
and the coordinates of the second electron and hole, 
and rh2- In our relative coordinate system these 
are represented by a vector to their centre of mass, 
i"2 = (fei +rhi — (re2 +rh2))/2, and a relative coordinate 
ra = re2 — r^2- The application of the dipole approxi- 
mation. Pop ~ requires the annihilating electron and 
hole to be at the same point in space. Thus, this relative 
coordinate is zero, and so 0be has the functional depen- 
dence 0BE(ri, r2, 0). For the bulk system, the centre of 
mass coordinates can be integrated over analytically. We 
then obtain a formula for the transition rates R into the 
bound states n, 



R71 — 



Pc 



(9) 



(iridr2(/)E,n(ri)*0BE(ri,r2,O)/(ri,r2) 



where we have introduced the function /(ri, r2) to model 
confinement in a quantum dot. We will discuss this func- 
tion in section III; for bulk, /(ri, r2) = 1. The additional 
factor of 2 in equation (jH]) comes from the fact that both 
excitons in the biexciton can recombine. 

The frequency w„ in equation ^ is the frequency 
associated with a specific transition. We neglect its 
momentum-dependence: The frequencies that contribute 
most to the emission spectrum will be close to the gap en- 
ergy Eg and hence we approximate a;„ — >■ Eg/h from here 
onwards. Our results confirm later that this treatment is 
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FIG. 1. Bulk emission spectra (overlap element as shown in 
Table HIl versus the energy difference of the surviving exciton 
with respect to its ground state, Ae/) for the two limiting 
mass ratios. The emitted photon energy corresponds to fku — 
Eg — B2X — IRy + Ae for irih S> rrie and = Eg — B2X — 
0.5Ry + Ae for ruh = nie, where B2X is the biexciton binding 
and a single exciton has energy Eg — IRy and Eg — 0.5Ry 
for these two mass ratios, respectively. Note that the peaks 
from the discrete bound states are delta functions and hence 
represented by arrows, labelled by their intensity. While it 
is clear that the continuous part of the spectrum decays to 
zero for large energies E, there is a possibility for a small but 
finite overlap for intermediate E. The number labelling the 
continuum is its integrated weight. 



indeed adequate. Thus the rate Rn{n) in equation (O is 
a function of quantum number n only. 

Since i?„ is the transition rate into a particular bound 
state, Rn and W{hhj) are related by 



(hio), (10) 



continuum 



where Wcontinuum(^) describes emission into continuum 
states. Wcontinuum(^) Can be obtained analagously to 
equation with the only difference that the sum over 
discrete quantum numbers n in equation ([5]) becomes an 
integral over the continuum states and the continuum ex- 
citon eigenfunctions are confluent hypergeometric func- 
tions. 



B. Results 

The emission spectra from bulk for the two limiting 
mass ratios, to^ — 1837me (corresponding to hydrogen), 
and me = nih are shown in Figure [1] As we have not 



TABLE II. Dimensionless overlap elements, 
2 J dridr2(^E,n(ri)*(^BE(ri, r2, 0), of excitonic bound (quan- 
tum number n) and continuum states with the biexciton 
wavefunction for various mass ratios. The mass ratios 
(rnh'-rrie) correspond to decreasing the rrih from 1837me (in 
hydrogen) by factors of 10. Individual elements are converged 
to approximately %0.5 accuracy. An estimate of the absolute 
accuracy including the uncertainty of the wavefunction 
cannot be given, but we tentatively estimate 10% accuracy 
from comparisons with a different variational form. 'Total' 
indicates the sum over the calculated bound and continuum 
states, and the result from equation ()12|) is shown in the last 
row. We note that for n > 2 the values drastically decrease 
compared to n = 1 (base change as shown in the first 
column); however, these higher n matrix elements depend 
very sensitively on the form of the wavefunction. 
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included broadening, the delta-function peaks for transi- 
tions into the bound states are indicated by arrows, and 
labelled with the corresponding delta- function coefficient. 
Figure [T] shows clearly that the transition into the ground 
state dominates. The spectra for the intermediate mass 
ratios show a similar behaviour. 

Individual overlap integrals for all different mass ratios 
are shown in Table [TTl which also contains the results ob- 
tained from equation (jl2l) . The discrepancy between the 
results from the individual overlaps ('total') and from 
equation (jl2p is a rough estimate of the numerical errors 
of the calculation. The agreement of these values shows 
that the results for the ground state are numerically sta- 
ble. Further error associated with the results are due to 
the approximate shape and the parameters involved in 
the wavefunction. A different variational form for hy- 
drogen gave bound state transition rates within 10% and 
continuum transition rates within about 30% of those in 
Table HT] While one can tell that the overlaps with the 
higher excited states decrease with n, the actual individ- 
ual numbers depend very sensitively on detailed prop- 
erties of the numerically determined wavefunction and 
hence should not be interpreted too closely. 

The general shape of the continuous part of our emis- 
sion spectrum should nevertheless be correct. The shape 
is notably different to the shape of the single exciton ab- 
sorption spectrum calculated by Elliott-. This single ex- 
citon absorption spectrum should correspond to a single 
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exciton emission spectrum, provided the single exciton 
emission is possibe: a single exciton in free space cannot 
recombine due to momentum conservation, but recombi- 
nation is possible when impurities are present and mo- 
mentum conservation is relaxed. Nevertheless, the con- 
tinuum part of the single exciton and biexciton emission 
spectra are different. Similar to single exciton absorp- 
tion, the biexciton continuum has a finite onset, but falls 
off rapidly with E, while the exciton absorption contin- 
uum line increases square-root-like for large energies. We 
expect the general trend of our spectrum to be correct, 
although we note that our continuum curve at energies 
of order E « IRy is numerically sensitive to details of 
the Monte Carlo wavefunction. In this energy range we 
sometimes obtain a small but finite overlap, which is due 
to the numerical integration being very sensitive to the 
variational form of the wavefunction in this energy range. 

The continuum emission decreases with E due to the 
decreasing wavefunction overlap of the initial state with 
the higher energy exciton state after photon emission. 
The fast decrease of the continuum emission in Figure 
[T]is especially notable for equal masses. Such a biexci- 
ton involving particles with equal masses is more loosely 
bound than one with a heavier hole. Loosely bound exci- 
tons have a smaller binding energy and thus less energy is 
available from the recombination. As the binding energy 
decreases the recoil momentum of the surviving exciton, 



K 



decreases towards zero, where recombination 



is prohibited by momentum conservation. The weaker 
fall-off of the continuum curve for the tighter bound biex- 
citon suggests that non-radiative processes, where the 
surviving exciton absorbs all the energy available from 
the annihilation, are more likely for heavier holes. When 
non-radiative decay is the main recombination mecha- 
nism, exciton lifetimes are dependent on charge carrier 
density or temperatureSiiii. 

The most notable feature of our calculation is the pre- 
dominance of the Is peak. This predominance implies 
that the surviving exciton is not affected by the recom- 
bination of the other-electron hole pair, even though we 
initially assume a strongly interacting bound state. Our 
model thus justifies treating excitons as non-interacting 
for the purpose of optical recombination. However, we 
speculate that interactions become more important if 
higher excited initial states are present. 

The lineshapes in emission spectra from biexcitons 
are often broadened due to scattering with phonons2Si21. 
This broadening smears out the different peaks from a 
variety of possible final states, and hence no direct com- 
parison with our spectra is possible. Similar problems 
prevent measurements of lifetimes in bulk^S. In addition, 
higher excited initial biexciton states may be present. 
However, our tentative spectra for bulk can be used as a 
reference point for other systems. 



III. LIFETIMES 

We now turn to the calculation of the biexcitonic life- 
times. We first discuss how to include polariton effects 
and present a general formula for the biexcitonic lifetime 
before turning to a brief discussion of lifetimes for a con- 
fined quantum dot. The results are also divided into a 
section on bulk and quantum dot systems. We compare 
to experimental and theoretical values in each of these 
subsections. 



A. Methods 

The overall rate for biexciton decay, which corresponds 
to the inverse biexcitonic lifetime, is made up from the 
individual transition rates in equation (|10p. 



1 
T 



d{huj)wihLu) 



Y.Rn + J dihuj)W, 



continuum 



(fie). (11) 



We note that the sum over excitonic states in equation 
([TT|) is complete and hence can be understood as a reso- 
lution of identity, and thus total biexciton decay rate can 
also be expressed in terms of the biexciton wavefunction 
only, 
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Pc 



2eo7r h?c^ 



2 y dr2dridr2 0BE(ri,r2,O)0BE(ri 



,0) 



(12) 



This expression should provide a check for the numerical 
errors associated with the integral. 

The photon generation rate in equation ^ does not 
yet include the polaritonic eigenstates. We mention ear- 
lier that these are important for nearly all semiconduc- 
tors, since the biexciton binding energy ^ Vic (the 
polariton splitting parameter) in almost all materials. 
The exact calculation of lifetime from a semiconductor 
would thus require the rate of polariton generation per 
unit energy 

k,S,(7 

(5(fiw - £;polariton(k, S,a-)), (13) 

where ^ g. {t) creates a polariton with momentum k 
and spin a in state s and -Epoiariton is the correspond- 
ing dispersion. The factor of one half avoids double- 
counting, since in the biexciton decay two polaritons are 
generated)^ In general the photon can be reabsorbed to 
form another biexciton and reemitted in many cycles. 
This would require a self-consistent treatment such as 
performed in Ref. 0, which we expect to be important 
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for materials with large Vlc/em- As we saw previously 
that the decay rate is dominated by the decay rate into 
the Is exciton, we assume that the created polariton is 
a superposition of a Is exciton and a photon and do not 
consider excited exciton states. We also use the resonant 
approximation for the polaritonic eigenstates, which we 
expect to cause an error of no more than 15%. With 
these approximations, our expression for the inverse life- 
time becomes 



1 = 2 

T 2 

[Eo{-Eo + 2Eg) + nl]^nl 



X |(*E,{n.}.Kl*BE>|' X {2(i?o - 2i?g)2 

} 



2{Eo-2EgmEo-2Eg)^ + n^^ 

73 1 



2 2 



TOq ireeQ Eg /h {hc)^ ' 



(14) 



Because lifetimes are longer and hence more exten- 
sively studied for confined systems, we wish to study 
the effect of confinement on the recombinative lifetime, 
and we consider a 3D quantum dot system. We assume 
a spherical quantum dot of radius d, where the photon 
wavelength A ^ d 3> aE,aB, where qe and as are rel- 
evant exciton and biexciton lengthscales. The exciton 
lengthscale as is the distance between the electron and 
hole in the surviving exciton, and the biexciton length- 
scale as is the distance between the centres of mass of 
the two excitons in the biexciton. Such a dot is large 
enough to only affect the centre of mass, but not the rela- 
tive wavefunction. The centre of mass wavefunctions can 
then be expressed in terms of spherical Bessel functions, 
and we obtain the constraint that the centre of mass 
wavefunction of the surviving exciton must also be in 
the ground state. Again, we use the Is exciton-polariton 
wavefunction; the corresponding lifetime is R(n = 1) in 
equation The effect of confinement is reflected in the 
function /(ri,r2) in equation ([9]), where 



/(ri,r2) 



4dsin(f^) 



1-2 



B. Results 

We have calculated lifetime estimates for a number of 
materials, and compare these to experimental and other 
theoretical lifetime estimates in this section. The pa- 
rameters required for the calculations of these lifetimes 
for CuCl, GaAs, ZnSe and InGaAs are given in table 
mil and table IIVI contains quantum dot values for those 
parameters which differ between bulk and quantum dot 
systems. Since the calculated lifetimes depend quite sen- 
sitively on the correct parameter and there is some vari- 
ation of parameters in the literature, this comparison is 
quite difficult. One example is e.g. the InGaAs mass ra- 
tio, where different groups use bulk or quantum well or 
interpolated quantum well mass ratios, or the experimen- 
tal GaAs lifetime. Here we quote an experimental value 



for GaAs bulk lifetime of 1.8 ns in table|V]merely for com- 
pleteness, since we expect the real lifetime to be around 
10-100 ps-". This estimate is based on the fact that the 
bulk biexciton lifetime should be smaller than the life- 
time in confined system like quantum dots or quantum 
wells, where the biexciton binding energy is higher^-'^ . For 
quantum wells, the biexcitonic lifetime is expected to be 
around half the excitonic lifetime^i^, which in Ref. ^ is 
quoted to be of order 200 ps, while Ref. [H contains an es- 
timate of 10 ps for the excitonic lifetime and Ivanov et al. 
0. quote lifetimes of order of magnitude 1-10 ps. There 
exists nevertheless a higher lifetime estimate of 1 ns^^ for 
quantum dots; however Ref [U cite potential impurities 
as a reason for lifetimes of this order of magnitude in 
quantum wells. Thus, the variety in experimental values 
makes the comparison sometimes somewhat difficult. 



Bulk 



We now discuss polariton lifetimes for bulk materials 
and the variation of the lifetime with mass ratio. 

Without including polariton effects the lifetime for 
CuCl according to equationlTil would be 20 and 100 ps for 
mass ratios 1.8 and 18, respectively. The experimental 
lifetimes are between 18 and 24p a'^^i^^ for a mass ratio 
mh/me ss 5, which shows that the inclusion of polari- 
tonic states (giving instead lifetimes 9ps and 47 ps for 
the same mass ratios) is indeed crucial. For the three 
semiconductors GaAs, ZnSe and CuCl experimental and 
our calculated lifetimes are shown in table |Vl We expect 
our method to work well for small ratios flc/em, while for 
large fic/sm a full bipolariton wavefunction and thus the 
inclusion of recursive exciton creation and annihilation 
processes is necessary. Additional errors in our estimate 
are associated with intrinsic errors of the wavefunction, 
the numerical integration, the uncertainty of the mate- 
rial parameters, the different geometry of the quantum 
dot and the fact that experimentally, the exciton might 
leak out. 

Indeed, our lifetime estimates for CuCl between 9 and 
47 ps for mass ratios 1.8 and 18, respectively, agree rea- 
sonably with the experimental value for the intermediate 
mass ratio of nih/me ~ 5. Ivanov et. ali^ also obtain ex- 
cellent agreement when using the bipolariton wavefunc- 
tion, but quote a poor lifetime estimate when using the 
giant oscillator model^ (4ps). The improved wavefunc- 
tion and a different weight in front of the process respon- 
sible for creating two lower state polaritons in equation 
[13] account for this difference. 

For ZnSe and GaAs the increased ratio ilc/sm indicates 
that polariton effects should matter more. Indeed, our 
lifetime estimates are lower than the experimental values 
by about factor of 10. 

Figure [5] shows plots of the dependence of the lifetime 
on mass ratio, with parameters for CuCl. The decrease 
of lifetimes on approaching the equal mass ratio limit in 
Figure HJa) is initially counterintuitive as electron and 
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TABLE III. Semiconductor parameters with references. Em is the biexciton binding energy. 



constant 




Ziiioe 




mo.6^£i'0.4 As 


Em/meV 


0.13^, 0.45^, 0.7± 0.2^" 


3.5*'^ 


33ii , 34i, 42^, 43^ 




Eg 


1 t-24 


2.8^ 


3.3^ 


0.7— 


e 


13.1^24 


8.1^ 


5.&L, 5.0^ 


14.3^, 14.03^ 


nih/mo 


0.52i 


1.7^ 


1.8ii, 2.0^ 


0.46^ 


me. /mo 


0.067— 


0.017^ 


0.5^'^ 


0.04^, 0.05^, 0.067^2 


Q.c/meV 


IS.&i 


100^ 


19li 


7 


Ep ^ 


25.7^ 


29.56 (from^) 


2.3 


21.7^ 



from bulk values, but inclusion of strain effects in Ifamiltonian 

Ep and Sic are related via the formula Oc = 2y^\/2npcv fi/ ('tioV i-KeLdofi^J and thus one can be derived from the other. The value 
of Ep in InGaAs where we do not quote a source was generated by this formula; note that the Qc of lOOfieV measured hyMS, is for a 
quantum dot in a nanocavity and thus not immediately applicable. For all other materials Ep is taken from the quoted source. 



TABLE IV. Quantum parameters with references. Em is the biexciton binding energy. 



constant 


CuCl 


Ino.6Gao.4 As 


Em/meV 


50aa (3nm), 51^, 60^ (3nm) 




mu/mo 


1.8 ^ 


0.125^^, 0.3^, 0.46^, 0.5^^ 



from bulk 

^ from interpolation between InAs and GaAs masses, using quantum well masses from™ and^ 
from bulk values, but inclusion of strain effects in Ifamiltonian 



TABLE V. Lifetime estimates for several semiconducting ma- 
terials and experimental values for comparison, rthcor corre- 
sponds to other theoretical calculations of biexciton lifetimes 
that are discussed in the text.. We expect the correct lifetime 
for GaAs to be of order 10-100 ps (see main text for discus- 
sion) . 



bulk 




quantum dot 


GaAs ZnSe 


CuCl 


CuCl Ino.6Gao.4 As 


r/ps 0.5-2.7 0.4-2.0 


8.8-47.3 


47.4-66.9 5.5-12.8 


Texp/pS 18002S 40^ 


18-27^1^ 


65^ 50C*' 


Tthcor /ps 


24''' 


500^2^ 



hole are much more likely to be at the same place for an- 
nihilation in the more tightly bound hydrogenic exciton. 
We find however that the probability for recombination 
is highest when the two excitons within the biexciton 
are far apart. This dependence on the mass ratio can 
be explained by the relative importance of the biexciton 
lengthscale qb (the distance between the centre of masses 
of the two excitons) and the exciton lengthscale ue (the 
distance between electron and hole within an exciton) in 
the overlap integral in equation The reason for the 

influence of the distance on the recombination rate 
comes from the overlap integral in equation the sin- 
gle exciton final state is delocalised. The overlap with 
the biexciton wavefunction is thus higher for a biexciton 
wavefunction in which the centres of mass of the exci- 
tons are more delocalised, i.e. where as is large. The 
lengthscale ratio as /as is inversely proportional to the 
ratio of biexcitonic to excitonic binding energies, and the 



400 
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1 



470 



460 - A 



450 



mil 3> me- 



5 10 15 20 25 30 

d (nm) 



10 100 1000 




10 15 20 25 30 

d (nm) 

(a) Lifetime versus mass ratio on (b) Lifetime versus dot size 
logscale 

FIG. 2. Left panel: Biexciton lifetimes (in ns) in 3D bulk, 
calculated with approximate values for a CuCl semiconductor, 
for different mass ratios as shown in Table [III Right panel: 
Biexciton lifetime for hole masses — 1837me and = m^ 
for different CuCl quantum dot radii. 



latter decreases towards the equal mass ratio limit for 
both two and three dimensional systems^iSS; thus, in 
the equal mass biexciton, the excitons are loosely bound 
and are further apart. 
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2. Quantum dot 

We now turn to lifetimes of confined systems, and dis- 
cuss both the impact of confinement on the lifetime as 
well as some calculated lifetimes for two materials. 

Figure [2Ib) shows the effect of confinement in a spher- 
ical quantum dot where confinement principally affects 
the centre of mass wavefunction. Stronger confinement 
increases the exciton localisation and thus the increase of 
lifetime with confinement shown in Figure [^fb) is consis- 
tent with the previous argument of increased annihilation 
probability for biexcitons made up of a loosely bound ex- 
citon pair. 

Unfortunately, experiments are not typically in the 
regime where confinement principally affects the centre 
of mass wavefunction, and so accurate numerical pre- 
dictions of the biexciton lifetimes in these cases require 
specific calculations for the particular structures. Such 
specific calculations have for example been performed 
for lno.6Gao.4As quantum dots. Using a configuration 
interaction method, Narvaez et al. [4^ obtain the ex- 
act experimental value of 0.5 ns for a 252x252x75 A 
Ino.6Gao.4As quantum dot, which agrees well with the 
experimental lifetime of 0.5 ns of a dot with height 2QA 
and a square base with length 150-200A'*^. Wimmer et 
al. [4^ perform a Monte Carlo optimisation of wavefunc- 
tions in anisotropic quantum dots for three materials, 
among them Ino.6Gao.4As, and obtain excellent lifetime 
estimates for all three materials. Since the exciton Bohr 
radius for this system is about 200A and the wavefunc- 
tion is thus strongly affected by confinement effects, our 
wavefunction is a inaccurate estimate for this system. 
Additionally, our calculation does not take into account 
strain effects and anisotropy of the hole masses, which 
are nonnegligible effects for InGaAs. For example, it is 
well known that these effects have a large impact on gyro- 
magnetic ratio^, where realistic k.p simulations specific 
to the material provide estimates adequate to the ex- 
perimental data^. Neglecting these effects and approx- 
imating this quantum dot by a spherical quantum dot 
of radius lOOA, we obtain a lifetime estimate of 6-13 ps, 
which is clearly inaccurate. 

For 3 nm CuCl quantum dots, the quantum dot radius 
is larger than the excitonic Bohr radius of 6. 8 A and thus 
our wavefunction should still be viable. Here, obtain a 
biexciton lifetime estimate of 47 -67 ps for our two limit- 
ing mass ratios, which is comparable to the experimental 
value of 65 ps. This shows that our method is indeed rea- 
sonable for the systems with small polariton effects, and 
can provide good estimates of the lifetime. 
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IV. CONCLUSIONS 

In this article we have investigated the lifetime and 
emission spectrum of a biexcitonic system in different 
mass limits. We employ the approach by Elliott in or- 
der to find the transition rates into excitonic states using 
Fermi's golden rule. We find that in our model, which 
assumes an initial biexciton in the ground state, the an- 
nihilation of one electron and hole is most likely to result 
in an exciton in its ground state. The predominance of 
this transition means that the two excitons which form 
the excitonic molecule can be treated as non-interacting 
for the purpose of optical recombination. Our estimates 
of lifetimes for different mass ratios are comparable to 
experimental values. 

In this calculation we have assumed a biexciton which 
was initially in its ground state. The presence of higher 
excited states could lead to a stronger overlap with higher 
excited final states, and thus to shorter lifetimes. We 
have also neglected the possibility of collision with other 
particles which could also lead to shorter lifetimes. We 
find that the biexciton lifetime slightly increases in con- 
fined quantum dots. Using an appropriately generated 
wavefunction, our method is also applicable for smaller 
quantum dots or quantum wells. 

Interesting extensions to fully model quantum dots 
could include strain effects and for example the impact 
of the wetting layer, for which experimental data is also 
available^. It would also be interesting to calculate life- 
times for pumped systems, where the timescale of the 
biexciton lifetimes is important for an accurate estima- 
tion of the efficiency of multi-exciton generation^SiZL. 

Our approach can be used as a benchmark for compar- 
ing approximative methods for systems at finites densi- 
ties. We also hope that the relative simplicity of the biex- 
citon emission spectrum as found here may imply that 
the development of further experimental and theoretical 
understanding of emission from the dense electron-hole 
plasma is within reach. 
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